Quadrupole deformation 7) of light A hypernuclei in constrained relativistic mean 
field model: shape evolution and shape polarization effect of A hyperon 
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The shapes of Hght normal nuclei and A hypernuclei are investigated in the (/3,7) deformation 
plane by using a newly developed constrained relativistic mean field (RMF) model. As examples, 
the results of some C, Mg, and Si nuclei are presented and discussed in details. We found that 
for normal nuclei the present RMF calculations and previous Skyrme-Hartree-Fock models predict 
similar trends of the shape evolution with the neutron number increasing. But some quantitative 
aspects from these two approaches, such as the depth of the minimum and the softness in the 7 
direction, differ a lot for several nuclei. For A hypernuclei, in most cases, the addition of a A hyperon 
alters slightly the location of the ground state minimum towards the direction of smaller /3 and softer 
7 in the potential energy surface E ~ (/3,7). There are three exceptions, namely, j^C, j^C, and 
i^Si in which the polarization effect of the additional A is so strong that the shapes of these three 
hypernuclei are drastically different from their corresponding core nuclei. 

PACS numbers: 21.60.Jz; 21.80.+a; 27.20.+n; 27.30.+t 
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I. INTRODUCTION 

Since the first observation of hypernuclei in 1953 [Q, 
a lot of experimental efforts have focused on the study 
of the spectroscopy of hypernuclei, see, for examples, 
Refs. for recent reviews. Due to the additional 

strangeness degree of freedom, a hyperon is free from 
nucleon's Pauli exclusion principle. Thus it can move 
deep inside the nuclei and may serve as an impurity for 
probing many nuclear properties that are not accessi- 
ble by normal methods. The study of hypernuclei can 
also provide detailed and accurate information about the 
hyperon-hyperon (YY) and hypcron-nuclcon (YN) inter- 
actions [ptoI which are important not only for the un- 
derstanding of hyper nuclear structure but also for the 
study of hyper matter and neutron stars . 

As an impurity in normal nuclei, a hyperon may induce 
many effects on the core nucleus, such as the shrinkage 
of the size ||lo|-[l^, the change of the shape which will 
be discussed later, the modification of its cluster struc- 
ture [l5[| , the occurrence of nucleon and hyperon skin 
or halo p^-p^, and the shift of neutron drip line to a 
neutron- rich side [0-|9| . 

The shape describes in an intuitive way the spatial 
density distribution of a quantum many-body system. 
Most of known nuclei are non-spherical and many are 
well-deformed as manifested by regular rotational spec- 
tra |2^. The shape-driven effect of valence nucleon(s) 
has been extensively studied in nuclear high-spin states, 
see, for examples, Refs. pl|-|2^. A well known example 
is that a nucleon occupying a high-j and low-f2 orbital 



would drive the nucleus to a more prolate shape. It is 
expected that the addition of a hyperon may also result 
in a shape polarization effect. Since the additional hy- 
peron is not restricted by the nucleons' Pauli exclusion 
principle, it tends to occupy the lowest s orbital, thus 
driving the core nucleus to be more spherical. If it oc- 
cupies a p orbital, a hyperon may enhance the nuclear 
deformation |25| . 

The self-consistent mean field models, either the non- 
relativistic [p6|"p9| or the relativistic ones [^0|-^, have 
been extensively used to reproduce the available hyper- 
nuclear data and/or make various predictions for hyper- 
nuclei. Up to now most of these studies focus on spherical 
systems. The first self-consistent mean field study of hy- 
pernuclei with an axially deformed Skyrme Hartree-Fock 
(SHF) model was finished by Zhou et al. |43 . It was 
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found that the core nuclei and the corresponding hyper- 
nuclei have similar deformations with the same sign 
which means that the shape polarization effect of the A 
hyperon is quite small. A further study within the same 
framework but with a microscopic AN force gives sim- 
ilar conclusions |Q. However, a relativistic mean field 
(RMF) study reveals that although in most cases the 
results are similar to the SHF calculations, there are 
indeed several exceptions, for examples, and ^^Si 
whose shapes change dramatically compared to their cor- 
responding core nuclei . The different results between 
the SHF and RMF calculations are attributed to the dif- 
ferent polarization effect of the additional A in these two 
approaches . In this sense, the experimental informa- 
tion related to shapes of hypernuclei is much desirable 
and would be used as a good benchmark for theories. 

The triaxiality is an important shape degree of free- 
dom in many nuclei. In an atomic nucleus with a sta- 
bly triaxial shape, the spontaneous broken chiral sym- 
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metry occurs p7|-p9[ and the wobbling motion is also ex- 
pected |l| . The triaxial deformation plays important 
roles in 7-soft nuclei and in nuclei in the transitional mass 
regions |p2|-3. The 7 deformation also changes consid- 
erably the local minima and the saddlepoint in the po- 
tential energy surface of heavy nuclei |5^ . Recently 
Win et al. have studied hypcrnuclci by using a SHF 
model with the triaxial degree of freedom included [ p7| . 
It is found that with an additional A hyperon no signif- 
icant change occurs for the nuclear shapes except that 
the potential energy surface becomes softer in the 7 di- 
rection. 

So far the relativistic description of hypernuclei is only 
restricted to spherical or axially deformed cases. It is the 
aim of this paper to investigate the triaxial deformations 
of the A hypernuclei and the shape polarization effect 
of the A hyperon in the RMF model, as a comparative 
study with the SHF calculations p^ . 

The paper is organized as follows. In Section ^ we 
briefly introduce the RMF model used in the hypernu- 
clear studies with emphasis on the extension of the RMF 
model to the inclusion of the hyperon and the triaxial- 
ity. In Section III, we present the calculated results for 
C, Mg, and Si isotopes and the corresponding hypernu- 
clei and discuss the shape polarization effect of the A 
hyperon. Finally a summary is given in Section [V. 



II. TRIAXIALLY DEFORMED RELATIVISTIC 
MEAN FIELD MODEL FOR HYPERNUCLEI 

In the relativistic mean field (RMF) model, the 
hadrons interact with each other via the exchange of cr, 
a;, p mesons and the photon. For hypernuclei, the RMF 
Lagrangian density can be written as: 



£ = £0 + -Ca, 



(1) 



where £0 is the standard RMF Lagrangian density de- 
scribing the nucleons and the couplings between nucleons 
and mesons [p8|-p2| and £a is that for the hyperon: 



(2) 



, /wAA , 
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where ttia is the mass of the A hyperon, g^rA and g^^A a-rc 
the coupling constants of the A hyperon with the scalar 
and vector meson fields, respectively. The last term rep- 
resents the tensor coupling between the A hyperon and 
the Lu field |6^. flf^i, is the field tensor of the lo field 
defined as fl^i, = — d^uj^. Couplings to the p me- 

son and the photon vanishes for A hypcrons which are 
neutral and isoscalar. 

Under the mean field approximation, the single particle 
Dirac equation for A hyperons reads: 

[d-p + fi (mA + Sa) + VX + Ta] ipM = eiiiAi, (3) 



with the scalar potential Sa = Qa-A'^-, the vector potential 
Va = .g^AW and the tensor potential: 



rp /wAA of- -\ 



(4) 



The potential energy surface (PES) is obtained by the 
constrained self-consistent calculation. 



E' ={H) + J2 7;Cn {{Qn 



n=l 



(5) 



where H is the RMF Hamiltonian, Qn's are the multi- 
pole operators to be constrained and is the dimension 
of the constraining space. The quadrupolc deformation 
parameters /3 and 7 arc calculated from the multipole 
moments of the baryon density distributions: 



(Q2 



-3(Q22) = 



7 = arctan 



V3(g: 



22/ 



{Q2 



(6) 



(7) 



where {Q20) and (Q22) are the quadrupole moments: 
{Q20) = I dr p{r){3z^-r^), 
{Q22) = f dr p{f){x^-y'). (8) 



For normal nuclei, the triaxially deformed RMF model 
has been developed based on expanding the nucleon 
Dirac spinor in a three-dimensional harmonic oscillator 
(3DII0) basis joj, |6^. As an alternative approach, in 
the present work, the RMF equations are solved in an 
axially deformed harmonic oscillator (ADHO) basis [ p6[ . 
We have modified the DIZ (or RMFAXIAL) code ||6^,]^ 
in order to allow the triaxial deformation and to include 
the hyperon. The basis wave functions are solutions of a 
Schrodinger equation with an ADHO potential: 

-^V^ + hi {cuy+u^lz')^ \a) = E^\a), (9) 



with r 



and 



\a) =<j>nAz)RZir)^c^p{imd)x. 
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(10) 



where a = {nz,nr,m, s) arc the asymptotic quantum 
numbers and Xs is for the spin. They are characterized 
by the basis deformation /3b- These basis states form an 
orthonormal complete set and can be used to expand any 
spinor wave functions irrespective with their symmetries. 

The projection of the total angular momentum on the 
symmetric z-axis K is not conserved due to the break- 
ing of the axial symmetry. The remaining symmetries 
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are discrete ones such as interchanges of the three axes. 
Namely, the system is invariant under the point group 

To describe the potentials and densities we use a dis- 
cretized three-dimensional mesh in the space. These 



mesh points are selected so that a Gaussian quadrature 
can be applied in the r and z directions, while an equally 
distributed mesh is used for the azimuthal angle. For 
convenience we use the Fourier expansion of the poten- 
tials V and densities p: 



2n 



00 ^ 

y]/n(2:,r)^C0s(2' 
n— 1 ^ 



nO), 



(11) 



with f = V or p. In an axially symmetric case only the n = term survives, the present code then returns back 
to DIZ. Most of the formulas for /„ are formally the same as the corresponding ones for /o, for example, the matrix 
element of the potential V between two basis states is: 



Vac,' = {n^_,nr,m, s\V\n'^,n'^,m' , s') 
1 



Ss.s' 
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V26K.K' Rn^nZK,n'^ + X! i^K'-K+2n,Q + fe'-/f-2n,o) -R^l™'"„^,„; 



n=l 



(12) 



where 



ris .rir ,n' 



dz 



rdr cpnA^)R'l{r)Vn{z,r)cp^,{z)Rl] (r), n = 0,1, 



(13) 



TABLE I. The calculated binding energies of ^®Si against A^f, 
the number of major shells for the Fermion basis. The pa- 
rameter set PKl is used. The deformation is constrained to 
P — 0.4 and 7 — 10°, 30°, and 50°, respectively. The basis for 
bosons are truncated up to Nb = 20. The unit for energies is 
MeV. 



Nf\ (/3,7) 


(0.4,10°) 


(0.4,30°) 


(0.4,50°) 


8 


-201.667 


-202.270 


-200.803 


10 


-201.523 


-202.093 


-200.633 


12 


-201.382 


-201.952 


-200.497 


14 


-201.325 


-201.907 


-200.468 


16 


-201.296 


-201.881 


-200.453 


18 


-201.294 


-201.878 


-200.453 



TABLE II. The RMF parameter sets used in the calculations. 





NN channel tua (MeV) 


Ra 


Rui 




PKl-Yl |j 


1 PKl jell 1115.6 


0.580 


0.620 


-1 


NLSH-A p 


5[ NLSIl|70|| 1115.6 


0.621 


0.667 


-1 



truncation errors arc less than 100 kcV for A^f > 12 and 
less than 30 keV for Np > 14 for these three points. 
Furthermore, when discussing the energy differences on 
which we focus in the present work, the truncation errors 
may even be less due to the cancellation. Therefore we 
use A^F = 14 in the following calculations. 



The Klein- Gordon equations for mesons are also solved 
by the basis expansion method, while the Coulomb field 
is solved by the Green's function method. 

In order to get a point on the PES with given deforma- 
tion parameters (/3,7), the axial deformation parameter 
of the basis is set to be, 



/3b = /3cos7. 



(14) 



The harmonic oscillator basis are truncated up to A^f 
fermion shells and A'b boson shells. The convergence of 
our method is checked for the nucleus ^^Si. Table | shows 
the calculated binding energies using different A^f- Be- 
cause the time consumption of the code is less affected 
by A'^B, we set A'^e = 20 which is big enough. Three typi- 
cal points on the PES are chosen, i.e., (/3,7) = (0.4,10°), 
(0.4,30°), (0.4,50°). In each column in Table | the same 
calculation is performed with A^f = 8, 10, ••• , 18. The 



III. RESULTS AND DISCUSSIONS 
A. Numerical details 

In this work we adopted for the Lagrangian density (|l|) 
two parameter sets which are listed in Table ||. For con- 
venience we give in Table || the coupling constants by 
three dimensionless quantities defined as R^ ~ 9a A /go-, 
Rui = guih/guj and i?„AA = fojAA/gujA- The effective 
interaction PKl-Yl is newly proposed by fitting to the 
experimental single-A binding energies and A spin-orbit 
splitting Based on the parameter set PKl |^ in 

the NN channel, PKl-Yl can reproduce the binding en- 
ergies of hypernuclei very well. For comparison we have 
also made calculations using the parameter set labeled as 
NLSH-A which is based on the NLSH parameter set for 
the nucleon-meson coupling constants [7^ and has been 
used in the axially symmetric RMF calculations for the 
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hypernuclei in Ref. p5| . 

Wc use a BCS scheme with constant gaps for the pair- 
ing. Following Ref. the pairing gaps are taken as: 

An = MeV, Ap = 4.8/^1/^ MeV. (15) 

In Ref. [ [72] , the fission barriers in actinides and super- 
heavy nuclei are calculated by using different pairing 
schemes and it is found that BCS calculations with con- 
stant pairing gaps do not provide an adequate descrip- 
tion of the fission barriers. The reason is if, e.g., a con- 
stant pairing strength G is used, the resulting pairing 
gap changes considerably with deformation because the 
density of single particle levels around the Fermi surface 
does so. This is a very important conclusion, especially 
for the study of fission barriers in heavy nuclei. For the 
shape evolution and shape polarization effect of the A hy- 
peron in light nuclei we investigate herejhe pairing does 
not play such a decisive role. In Ref. results from a 
constant pairing gap and a constant pairing strength G 
are compared for the light hypernuclei in the RMF-I-BCS 
model. It is shown that most of the results are only 
slightly changed and the main features and conclusions 
remain the same. However, if the potential energy surface 
is rather soft, the minimum or minima in the potential 
energy surface would be different when different treat- 
ments of the pairing correlation are used. In such cases 
the configuration mixing effects must be also included, 
as discussed later. 

The center of mass correction is included either phe- 
nomenologically (for the parameter sets NLSH and 
NLSH-A) or microscopically (for the parameter sets PKl 
and PKl-Yl). Note that neither the tensor force nor 
the center of mass correction has significant influence on 
shapes of the nuclei investigated in this work. For normal 
nuclei and single-A hypernuclei, they only shift the PES's 
by roughly a few MeV as a whole. But when discussing 
the absolute value of the energies, their contributions are 
certainly not negligible. 

Using this triaxially deformed RMF-I-BCS method, we 
have calculated the PES's of even-even C, O, Ne, Mg, Si 
and S isotopes. Then by adding one A hyperon, we have 
investigated in details the changes of the PES's of the cor- 
responding hypernuclei. Next we take carbon and silicon 
isotopes as examples and examine the shape evolution 
in these two isotopic chains and the shape polarization 
effect of A hyperon. Some results of ^^Mg and ^^Mg are 
also given for making a comparison with ^^Si and a'SI. 



B. Carbon hypernuclei 

The ground state properties of carbon isotopes have 
been studied extensively with the axially deformed 
RMF [||-^ and SHE models |7|l and the triaxially de- 
formed SHE model |7^. Here we generalize the RMF 
calculations by considering two additional degrees of free- 
dom: the triaxial degree of freedom for nuclear shape and 



the hyperon degree of freedom. We have performed con- 
strained RMF calculations in the (/3, 7) plane for even- 
even carbon isotopes from ^''C to the drip line nucleus 
^^C as well as for the hyper counterparts. Note that the 
study of carbon isotopes with triaxially deformed SHE 
models for normal nuclei fr\ and hypernuclei js^ ] are 
both available, while that with the RMF model is still 
absent. 

The calculated ground state deformation parameters 
are summarized in Table HI , together with the root mean 
square radii, binding energies and single A separation 
energies. The single A separation energy is defined as 
the energy difference between a hypernucleus and the 
corresponding core nucleus: 



B 



(A+l 



A U 



C) = ^ (^C) - E {^+^G) 



(16) 



Firstly, the contraction due to the additional A are ob- 
served for all the nuclei studies here, which is a manifes- 
tation of the glue like effect of the hyperon iQ . Secondly, 
the radii of ^"C and are much larger than those of 
the neighbouring nuclei with two more neutrons. This is 
due to that they are close to the proton drip line. We 
mention that the earlier calculations with an axially de- 
formed RMF-f-BCS model also show a similar trend [ [73| . 
Finally, because the A hyperon always occupies the low- 
est orbital generated by the mean field, the single A sep- 
aration energy i?A can be seen as a measure of the depth 
and shape of the potential felt by the hyperon. In gen- 
eral this quantity should increase as the nucleon number 
increasing because the dept h of the potential increases. 



This is clearly seen in Table III, 

In Figure ^ we present the calculated PES's of car- 
bon isotopes, together with those of the corresponding 
carbon hypernuclei with one additional A hyperon. The 
locations of the ground states are denoted by open trian- 
gles. To unify the energy scales, we only show the relative 
energies with respect to the ground state. The contours 
join the points on the PES with the same energy. The 
energy difference between two neighbouring contours is 
0.15 MeV. 

Many nuclei show two minima in both the prolate and 
the oblate sides of the PES's in axially deformed calcula- 
tions. In the triaxially deformed calculations, at least one 
of them becomes a saddle point. In such cases, the rel- 
ative energy differences among the minimum, the saddle 
point and the spherical configuration are key quantities 
characterizing the PES. We list in Table IV the ener- 



gies of the minimum or the saddle point in the prolate 
^'prolate or oblate sides -Eobiatc as well as that of the spher- 
ical configuration i?sphoricai, with respect to the ground 
state energy. The results calculated with PKl-Yl and 
NLSH-A parameter sets are both presented. 

It is convenient for further discussions to define the fol- 
lowing two quantities. One is the deformation energy de- 
fined as the energy difference between the spherical shape 
and the ground state, i.e., -Bdef = |£^sphericai - £^ground|- 
i?dct characterizes the driving force to deformation quali- 
tatively. The other is the prolate-oblate energy difference 
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TABLE III. The deformation parameters P and 7, root mean square radii r and binding energies E of the carbon nuclei and 
hypernuclei calculated with parameter set PKl-Yl. The subscripts n, p, A and tot represent the corresponding quantities for 
neutron, proton, A hyperon, and the whole nucleus, respectively. The single-A separation energies BAcai are also presented for 
the hypernuclei. 



Nucleus 






P 










7 (deg) 




r. m. 


s. radii (fm) 


Energies 


(MeV) 








/3p 




/3a 


Plot 


7n 


7p 


7a 


7tot 


rn 




r-A 


rtot 




BAcal 




0.340 





255 









284 













2.53 


2 


87 






2.74 


56.404 






0.260 





216 





060 





218 














2.47 


2 


79 


2 


43 


2.65 


65.652 


9.250 


12^ 


0.207 





213 









210 


60.0 


60.0 




60.0 


2.53 


2 


57 






2.55 


87.764 




fc 
































2.46 


2 


50 


2 


27 


2.47 


99.752 


11.990 


"c 






























2.74 


2 


52 






2.65 


104.943 




































2.73 


2 


50 


2 


37 


2.62 


117.152 


12.209 




0.225 





179 









211 













2.90 


2 


45 






2.74 


110.398 






0.202 





158 





068 





183 














2.88 


2 


43 


2 


40 


2.70 


123.422 


13.024 


18c 


0.258 





272 









261 


35.0 


40.6 




36.4 


3.23 


2 


57 






3.03 


114.874 




fc 


0.246 





256 





106 





243 


34.8 


40.5 


35.7 


36.2 


3.20 


2 


54 


2 


44 


2.97 


128.434 


13.560 


20 c 


0.275 





308 









282 


60.0 


60.0 




60.0 


3.47 


2 


63 






3.24 


119.601 




A 


0.259 





291 





119 





261 


60.0 


60.0 


60.0 


60.0 


3.43 


2 


60 


2 


47 


3.17 


133.591 


13.990 


22 c 


0.178 





268 









193 


60.0 


60.0 




60.0 


3.64 


2 


65 






3.40 


120.735 




fc 
































3.55 


2 


58 


2 


44 


3.28 


135.340 


14.605 
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FIG. 1. (Color online) The potential energy surfaces of carbon hypernuclei and the corresponding core nuclei in the (/3,7) 
plane, calculated using the PKl-Yl parameter set. The energies are normalized with respect to the binding energy of the 
absolute minimum. The contours join the points with the same energy. The contour interval is 0.15 MeV. The ground states 
are denoted by open triangles. 
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TABLE IV. The energies of the minimum or the saddle point 
in the prolate or oblate sides, and the energy of the spher- 
ical configuration, with respect to the ground state for car- 
bon nuclei and hypernuclei, calculated with the PKl-Yl and 
NLSH-A parameter sets. The energies are in MeV. 







PKl-Yl 






NLSH-A 


-^'prolate 


-Boblatc E 


sphcrcial 


-^'prolate 


-Eoblatc 


-S'splicrcial 







0.74 


0.95 





0.53 


0.70 







0.34 


0.37 





0.25 


0.29 


12q 


- 





0.26 


- 





0.11 


aC 
















































16c 





0.18 


0.28 





0.20 


0.34 







0.12 


0.18 





0.14 


0.22 


18c 


0.30 


0.04 


1.95 


0.22 


0.04 


1.86 


fc 


0.14 


0.05 


1.57 


0.08 


0.14 


1.52 


20c 


2.28 





2.88 


2.15 





2.75 


21^1 
A 


1.83 





2.28 


1.74 





2.14 


22c 







0.54 







0.25 


fc 

















defined as the energy difference between the minimum 
in tlie prolate (oblate) side and the saddle point in the 
oblate (prolate) side for an axially deformed nucleus or 
the energy difference between two saddle points in the 
oblate and prolate sides for a triaxially deformed nucleus, 
i.e., Epo = l^^proiatG - -Eobiatol- Apparently it is only weh 
defined if at least one saddle point exists in the prolate 
or oblate sides of the PES. Epo characterizes the softness 
of the PES in the 7 direction in most cases. That is, the 
larger Epo is, the steeper the PES in the 7 direction is. 



1. Shape evolution of carbon isotopes 

First let us examine the shape evolution of the car- 
bon isotopes. The isotopic dependence of the deforma- 
tion with triaxiality has been investigated with the SHE 
model in Ref. 0, where the PES's of even-even carbon 
nuclei in the (/3, 7) deformation plane are presented and 
discussed. Next one can find that there are some new 
features from the RME calculations. 

The energy minimum of -'^^C situated at /3 w 0.28 and is 
rather stable against the triaxial distortion. The energy 
difference between the ground state and the saddle point 
on the oblate side Epo is 0.74 MeV, accounts for more 
than 1% of the total binding energy. However, this ciuan- 
tity is less than 40 keV in the SHE calculations |57[ 0, 
indicating that is rather 7-soft from the SHE calcu- 
lations. To examine the parameter dependence of our re- 
sults, we also performed the RME calculation using the 
NLSH parameter set. The resulting Epo is 0.53 MeV. 
Results for from other RME parameter sets can also 
be found in Ref. ||7^ where Epo is always larger than 
0.5 MeV. Thus we conclude that the PES of is ap- 
parently softer in SHE than in RME in the 7 degree of 



freedom. 

With two additional neutrons in the sd shell, ^^C is 
driven to be oblate deformed. The driving force is so 
weak that the deformation energy E^ef is only 0.26 MeV 
and 0.11 MeV from the PKl and NLSH parameter sets, 
respectively. The SHE calculation gives a E^ef as large as 
1 MeV 0. This nucleus is empirically known to be 
oblate in its ground state from the inelastic scattering ex- 
periments [[78|-p0|, which is in consistent with our result. 
It is appropriate to compare the depth of the energy min- 
imum obtained from different mean field models. Most 
of the recent SHE calculations of the carbon isotopes use 
a recipe that the spin-orbit interaction is reduced to 60% 
of its original strength in order to reproduce the oblate 
shape of ^^C. This prescription is rather arbitrary and 
it is known that the depth of the energy minimum with 
respect to the spherical configuration is sensitive to the 
spin-orbit interaction (see Eig. 1 in ReL [|7^). In ReL 
it is shown that using this reduction factor as an ad- 
justable parameter, different shapes of a A hypernucleus 
could be predicted. In contrast with the SHE method, 
the deformation of ^^C is correctly reproduced without 
any adjustment of parameters in the RME model here 
and in Ref. ||7^ . It seems that the reduction of the spin- 
orbit interaction introduced in these SHE calculations is 
a little too strong, if the RME results are reliable. This 
needs to be further explored. 

With the neutron number = 8 which is magic, ^^C 
is predicted to be spherical. When two and four more 
neutrons are added to the sd shell, ^^C and turn to 
be prolate and triaxially deformed, respectively. Eor ^^C, 
-Epo is 0.18 MeV with the PKl parameter set and 0.20 
MeV with the NLSH parameter set while it is predicted 
to be 0.61 MeV in the SHE calculation Contrary 
to the result of ^"C, the energy minimum of ^^C is much 
deeper in SHE than in RME. The softness of the PES of 
^^C is also obtained by using three different mean field 
models in ReL ^^C is the only nucleus with a tri- 

axial deformation in carbon isotopes. There are saddle 
points at both the prolate and the oblate sides. The 
RME calculation predicts that the oblate one is lower 
than the prolate one by 0.26 MeV, while from SHE calcu- 
lations the prolate one is lower (cL Eig. ^Ihere and Eig. 5 
in ReL @). The PES's of ^^C and ^^C are rather 7- 
soft. The ground state deformations of such nuclei with 
extremely soft PES's may not be well described in the 
mean field level because the ground state wave functions 
are always a strongly correlated superposition of different 
shapes with nearly the same energies. This suggests that 
further investigations of these nuclei should include be- 
yond mean field effects by using the generator coordinate 
method | |8^ , |83[ |. 

The shape evolves again to be oblate for ^"C and ^^C 
with the latter to be the last bound nucleus within the 
neutron drip line of carbon isotopes. ^"C is strongly de- 
formed and its PES is the steepest one among the carbon 
isotopes investigated here in either /? or 7 directions. ^^C 
is suggested to be a halo nucleus according to the mea- 
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surcd large enhancement of the reaction cross section for 
it compared to those for neighboring carbon isotopes Q . 
We certainly can not reproduce the halo structure for ^^C 
because in the present work a Harmonic Oscillator basis 
is used which is not able to giv e the large spatial density 
distributions in halo nuclei p5|-p8[. 

The shape evolution of normal carbon nuclei can be 
roughly explained by examining the shell structure. One 
can see from the Nilsson diagram that 10 and 14 are pro- 
late and oblate magic numbers respectively |^^, which 
are responsible for the prolate shape of ^^C and oblate 
shape of 20c. 

In the above discussions we see that the RMF and SHF 
models predict the same trend of the shape evolution 
with the neutron number increasing. But some quanti- 
tative aspects from these two approaches, such as the 
depth of the minimum and the softness in the 7 degree 
of freedom, differ a lot for some nuclei. 



2. Shape polarization effect of A hyperon in carbon 
hypernuclei 

Next let us discuss the shape evolution in carbon hy- 
pernuclei and the shape polarization effect of the A hy- 
peron. As is indicated in Ref. the influence of a 
hyperon on the PES can be as large as 1 MeV, the shape 
may evolve in a different way compared with normal ones. 

Roughly speaking, the deformations of the carbon hy- 
pernuclei are similar to their corresponding core nuclei, 
with two exceptions and ^"^C. These two nuclei be- 
come spherical while the corresponding core nuclei ^^C 
and 22 C are both oblate. The spherical shape of ^'^C 
has been predicted in an axially deformed RMF cal- 
culation [Q. In our triaxially deformed RMF model 
this is confirmed in the 7) deformation plane. As 
we discussed earlier, the deformation energy of ^^C is 
quite small. From Fig. |l| one can see that at the oblate 
(7 = 60°) edge ^^C is very soft in the (3 direction. The 
additional A drives to be spherical. ^'^C becomes 
spherical but it is also very soft in the /3 direction at the 
oblate edge. The calculation with NLSH-A parameter 
set also predicts spherical ground states for ^'^C and ^'^C. 

In the SHF calculations, the PES of with a prolate 
shape is so soft in the 7 degree of freedom that one addi- 
tional A drives the shape of to be oblate In the 
present work, is still prolate with a smaller /3 com- 
pared to ^°C. The addition of a A hyperon only makes 
the PES of a little softer than that of i°C. Similar 
situation holds for ]^C. For the only triaxially deformed 
carbon nucleus, ^'^C, it is observed that the A hyperon 
also makes the PES of j^C softer and the 7 softness in- 
creases more towards the prolate direction. Since ^^C 
is strongly oblate deformed, the addition of a A hyperon 
does not change its shape much. Interestingly, with one A 
added, the PES of ^^C becomes stiffer around the spher- 
ical minimum compared to its core nucleus ^*C. This is 
also due to the spherical-driven effect of the A hyperon. 



TABLE V. The energies of the minimum or the saddle point 
in the prolate or oblate sides, and the energy of the spher- 
ical configuration, with respect to the ground state for sih- 
con nuclei and hypernuclei, calculated with the PKl-Yl and 
NLSH-A parameter sets. The unit of the energies is MeV. 
The results for ^^Mg and i'^Mg are also presented. 



PKl-Yl NLSH-A 





-Eoblatc 


-S^prolatc 


-^'spherical 


-Eoblatc 


-E'prolatc 


-^'spherical 


''Si 


- 


- 





- 


- 





23 

A til 


- 


- 





- 


- 





'''Si 


0.63 





1.17 


0.17 





0.62 


A Si 


0.40 





0.77 


0.07 





0.33 


'«Si 


1.02 





2.20 


0.45 





1.32 


fSi 


0.79 





1.55 


0.41 





0.87 


'»Si 







1.04 







0.40 


A 







0.10 














0.17 


0.53 





0.03 


0.32 


31c- 

A -SI 


<0.01 





0.12 


0.05 





0.14 


3'Si 





0.61 


0.66 





0.49 


0.52 


33c; 





0.41 


0.44 





0.32 


0.34 


3*Si 
















35 c; 
A t*! 
















'«Mg 


0.34 


0.04 


1.38 


0.32 


0.01 


1.08 


fMg 


0.29 





0.89 


0.33 





0.70 



The softness of the PES of a nucleus in the 7 direction 
can be measured by Epo. For example, it is 0.74 MeV for 
and 0.34 MeV for ]^ C, which means that the PES of 
the latter is much softer than the former. However, this 
difference for ]!JC is only 0.06 MeV smaller than that for 
the corresponding core nucleus, which indicates a very 
tender change. 

From the above discussion it is seen that the shape 
evolution of carbon isotopes is modified due to the addi- 
tional A hyperon. On one hand, the transition point from 
the deformed shape to the spherical shape is shifted from 
^"'C in normal nuclei to in hypernuclei. On the other 
hand, an abrupt change from a strongly oblate shape to 
a spherical shape is observed for ^■^C. Thus the sponta- 
neous symmetry breaking effect in carbon hypernuclei is 
very different from that in normal nuclei. 



C. Silicon hypernuclei 

As another example for the study of the shape evolu- 
tion between two shell closures, we present in Fig. ^ the 
PES's of silicon isotopes and the corresponding one-A 
hypernuclei in the (/?, 7) plane. Because in silicon iso- 
topes the proton number is 14 which energetically favors 
the oblate shape, the competition between the neutron 
and proton deformation driving forces may produce var- 
ious types of PES's. When a A hyperon is added, subtle 
changes are expected. 

We list the energies of the minimum or the saddle point 
at the prolate and oblate sides and the spherical config- 
uration, with respect to the ground state in Table M. As 



8 




0.3 0.4 0.5 



1.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 0.5 

PCOSy 



FIG. 2. (Color online) The potential energy surfaces of silicon hypernuclei and the corresponding core nuclei in the 7) plane, 
calculated using the PKl-Yl parameter set. The energies are normalized with respect to the binding energy of the absolute 
minimum. The contours join the points with the same energy. The contour interval is 0.3 MeV. The ground states are denoted 
by open triangles. 




0.1 0.2 0.3 0.4 0.5 0.1 
pcoSy 



0.3 0.4 



FIG. 3. (Color online) The potential energy surfaces of Mg 
and A'Mg in the (/3, 7) plane calculated using the PKl-Yl 
parameter set. The energies are normalized with respect to 
the binding energy of the absolute minimum. The contours 
join the points with the same energy. The contour interval is 
0.3 MeV. The ground states are denoted by open triangles. 



the muTor nucleus of ^^Si, the results for ^^Mg are also 
presented. For silicon isotopes the parameter dependence 
of the results is a bit larger than that for carbon, but the 
results from these two different parameter sets are still 
in consistent with each other. 

Starting from the spherical magic number nucleus ^^Si, 



the shape evolves to be prolate in ^^Si and ^^Si. These 
two nuclei are rather 7 soft, similar as ^^C and ^'^C. In the 
previous SHF calculations |57 an oblate ground state was 
obtained for ^^Si. Following Ref. |^ we also calculate 
26Mg which is the mirror nucleus of ^^Si. The PES of 
^^Mg is presented in Fig. |[ The ground state of ^^Mg 
is triaxially deformed. Although the calculated PES's 
for ^^Mg from the RMF and SHF calculations are both 
flat against the 7 deformation, our RMF model predicts 
a lower saddle point at the prolate side while the SHF 
method does at the oblate side (cf. Fig.|3|hcre and Fig. 12 
in Ref. [0). For hypernuclei ^^Si and ^''Si, the values 
of the deformation parameter j3 are a bit smaller than 
those of the corresponding core nuclei, respectively and 
the PES's become softer along the 7 direction. 

With two more neutrons in the sd shell, ^*Si turns 
to be oblate with £^dcf as large as 1.04 MeV. Adding 
one A hyperon results in a shape coexistence. Namely, 
the ground state is still oblate, but the energy of the 
spherical configuration with respect to the ground state 
is lowered to as small as 0.1 MeV. The barrier height be- 
tween the two local minima is only about 0.25 MeV. The 
results discussed above are obtained from the parameter 
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set PKl. Note that in the axially deformed RMF calcula- 
tion in Ref. [Q, ^^Si is predicted to be spherical with the 
NLSH-A parameter set. When the NLSH-A parameter 
set is used in our calculation, we obtain the same conclu- 
sion as Ref. [|5|. As is seen in Table 0, for ^*Si E'dcf is 
0.40 MeV, but for ^^Si E^d vanishes which means that 
^^Si is spherical. Therefore the prediction of the shape 
of ^^Si is a bit parameter dependent. 

The PES of "^"Si is almost the softest one among the 
nuclei investigated in the present work. Including an ad- 
ditional A hyperon, the energy of ^^Si is almost irrelevant 
with the deformation if /3 < 0.2. The ground state moves 
from the oblate side to the prolate side. But the shift of 
the minimum does not mean much because the coherent 
superposition of the different shapes should be considered 
for the actual ground state. ^^Si is also oblate and the 
PES is soft near the ground state. The A hyperon softens 
slightly the PES of ^'^Si. Filling completely the neutron 
sd shell, stable spherical shapes are again obtained for 
34Si and f Si. 



IV. SUMMARY 

We developed a triaxially deformed RMF model for 
hypernuclei. Different from previous RMF calculations 
for normal nuclei, in the present work, the RMF equa- 
tions are solved in an axially deformed harmonic oscil- 
lator (ADHO) basis. The convergence of the calculated 
results against the basis truncation is studied and it is 
shown that a reasonably large ADHO basis is able to pro- 
vide desired accuracy in the triaxial RMF calculations. 

The shapes of C, O, Ne, Mg, Si and S A hypernuclei are 
investigated in the (/3, 7) deformation plane by using this 
newly developed constrained RMF-I-BCS model with the 
parameter sets PKl-Yl and NLSH-A. As examples, the 
results of some C, Mg, and Si A hypernuclei are presented 
and we discussed in details the shape evolution of light 
normal nuclei and hypernuclei and the shape polarization 
effect of the A hyperon. 

It is found that for normal nuclei the present RMF 
model and previous Skyrme-Hartree-Fock models predict 
similar trends of the shape evolution with the neutron 
number increasing. But some quantitative aspects from 
these two approaches, such as the depth of the minimum 
and 7 softness differ a lot for several nuclei. 

For A hypernuclei, in most cases, the addition of a A 
hyperon alters slightly the location of the ground state 
minimum towards the direction of smaller /? and softer 7 
in the potential energy surface (PES) E ^ (/?,7)- There 
exist two exceptions in carbon isotopes, namely, and 



^^C in which the polarization effect of the additional A 
is so strong that it drives these nuclei from the oblate 
shape to spheres. Shape changes also occur in silicon iso- 
topes. Ahhough ^^Si is oblate from both the PKl-Yl 
and NLSH-A parameter sets, the prediction of the shape 
of ^^Si is parameter dependent. ^^Si is spherical with 
the NLSH-A parameter set but it is still oblate with the 
PKl-Yl parameter set. Compared to the core nucleus, 
shape change also happens in ^^Si, from an oblate shape 
in 30 Si to a prolate one in fSi. But the PES's of ^^Si and 
^^Si are rather soft. The ground state deformation of a 
nucleus with such an extremely soft PES may not be well 
described in the mean field level because the ground state 
wave function should be a strongly correlated superpo- 
sition of different shapes with nearly the same energies. 
This suggests that further investigations of these nuclei 
should include beyond mean field effects by using, for 
example, the generator coordinate method. 

Finally we note that since different predictions about 
the shape polarization effect of the A hyperon are made 
by different models and some times even by different ef- 
fective interactions within the same model, the exper- 
imental information related to shapes of hypernuclei is 
highly desired and would be used as a good benchmark 
for theoretical models. In Ref. |90 , Yao et al. studied the 



impurity effect of the A hyperon on collective excitations 
of nuclei based on potential energy surfaces calculated 
from the Skyrme HF model. It is found that the A hy- 
peron stretches the ground state band of the core of ^^Mg 
and reduces the B{E2 : 2^ 0^) value considerably due 
to a softening effect of A on the potential energy surface. 
We would expect more profound effects in the collective 
spectrum from an additional A if it changes even the 
shape of a nucleus, e.g., from a prolate shape in ^^C to 
a spherical one in ^^'^C. We expect that with the new or 
updated experimental facilities in J-PARC or JLab such 
measurements may become possible in the near future. 
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